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A recently discovered feature of hadron production in nuclear collisions is 
the large stopping of the participating nucleons in hadron-nucleus and nucleus- 
nucleus collisions. Experimental data demonstrating this effect have been pre- 
sented in Jl|,| . 

Multistring fragmentation models like the Dual Parton Model (DPM) or sim- 
ilar models did originally not show this enhanced stopping in nuclear collisions. 
Therefore, in order to incorporate the effect into multistring fragmentation mod- 
els new diquark breaking DPM-diagrams acting in hadron-nucleus and nucleus- 
nucleus collisions were proposed by Kharzeev || and Capella and Kopeliovich 
[fi} and investigated in detail by Capella and collaborators (|||. Similar ideas 
were discussed by Vance and GyulassyJrJ and by Casado ||. 

The Monte Carlo implementation into DPMJET-II.5 of the new diquark 
breaking diagrams of Kharzeev ||] and Capella and Kopeliovich [0 was first 
discussed by Ranft || . The implementation into DPMJET-III (l^] of these 
diagrams differs somewhat from || and is described here. 

There are two possibilities for the first fragmentation step of a diquark. Either 
we get in the first step a baryon, which contains both quarks of the diquark and 
the string junction or in the first step a meson is produced containing only one of 
the two quarks and the baryon is produced in one of the following fragmentation 
steps. This mechanism was implemented under the name popcorn fragmentation 
in the Lund chain fragmentation model JETSET |n||^] which is presently used 
in DPMJET. The popcorn mechanism alone is not enough to explain the baryon 
stopping observed experimentally in hadron-nucleus and nucleus-nucleus colli- 
sions 

In Ref.|D we describe these new diquark breaking DPM-diagrams in detail. 
This will not be repeated here. The two important diagrams are 

(i) GSQBS, the Glauber sea quark mechanism of baryon stopping, this dia- 
gram acts in nuclear collisions already at low energy. 

(ii) USQBS, the unitary sea quark mechanism of baryon stopping, this mech- 
anism leads to baryon stopping also in proton-proton collisions at collider and 
cosmic ray energies. 
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In DPMJET-III we first construct the system of parton chains according 
to the model without the diquark breaking diagrams. Having this we search 
for situations (i) as plotted on the left hand side of Fig. |l| or situations (ii) 
(not plotted), where the left lower diquark is replaced by an antiquark. On 
the left hand side in Fig. |l| we have a diquark-quark chain and a seaquark- 
diquark chain with the (upper) diquark and seaquark belonging to the same 
projectile hadron and also the (lower) diquark and valence-quark belonging to 
the same target-nucleon. In situation (ii) we have a diquark-quark chain and a 
seaquark-anti-seaquark chain with the (upper) diquark and seaquark belonging 
to the same projectile hadron and also the (lower) anti-seaquark and valence- 
quark belonging to the same target-nucleon. The chain system is transformed 
as plotted on the right hand side of Fig.|l|. The projectile diquark is split and the 
two resulting quarks come to the upper ends of both chains and the projectile 
seaquark goes into the middle of the second chain and determines the position 
where the baryon is produced. The sea quarks in Fig. [j] might be Glauber sea 
quarks or unitary sea quarks. 




Fig. 1. Diquark breaking in an original diquark-quark and seaquark-diquark chain 
system. 

Besides the situations already discussed we consider also the ones where 
projectile and target are exchanged as well as situations with anti-diquarks. 

We split the diquark by sampling for one of the two resulting valence quarks 
(randomly choosen) a normal valence quark x v \ and give to the second valence 
quark x V 2 — Xd~ x v \. For each of the new diquark breaking diagrams (GSQBS 
and USQBS) we have to introduce a new parameter. These parameters give the 
probability for the diquark breaking mechanisms to occur, given a suitable sea 
quark is available and given that the diquark breaking mechanism is kinemati- 
cally allowed. For an original diquark-quark chain of small invariant mass, which 
originally just fragments into two hadrons, the diquark breaking is often not al- 
lowed. The optimum values of the new parameters are determined by comparing 
DPMJET-III with experimental data on net-baryon distributions. We obtain 
for the GSQBS parameters the value 0.6 and use the same value for the USQBS 
parameter. 

Introducing the new baryon stopping mechanisms into DPMJET we get an 
significant modification of the model in different sectors: (i)The Feynman-x dis- 
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tributions of leading protons in proton-proton and proton-nucleus collisions. 
The leading particle production is especially important for the cosmic ray cas- 
cade simulation. (ii)The net-p (p - p) and net-A (A - A) rapidity distributions 
in hadron-nucleus and nucleus-nucleus collisions. These are the data on the 
enhanced baryon stopping mentioned already above. (iii)The production of hy- 
perons and anti-hyperons in nuclear collisions. We present here examples for (i) 
and (ii). 

In Fig. H we compare the distribution in the energy fraction xi a b carried by 
the leading proton. The data are photoproduction and DIS measurements from 



the HERA collider at y/s « 200 GeV g3). We compare to DPMJET-III for p-p 



collisions at y\s = 200 GeV. The forward production of leading protons is not 
expected to depend strongly on the reaction channel. We present the DPMJET- 
III distributions for the models with and without the diquark breaking diagrams 
(No Sto in the plot). The new diagrams modify the distributions mainly at 
intermediate x^ab values where, unfortunately, no experimental data exist at 
present. Therefore, further conclusions cannot be drawn at this point. 
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Fig. 2. Energy fraction xi a t carried by the leading proton. The data are photopro- 
duction and DIS measurements at %/s ~ 200 GeV compared to DPMJET-III with 
and without (No Sto) the diquark breaking diagrams for p-p collisions at y/s = 200 
GeV. 



In Fig. |3| we compare the net-proton distributions according to the models 
with and without the diquark breaking diagrams with data in p-Au collisions 
[|J. The dip at central rapidity, which occurs in the model without the baryon 
stopping diagrams is filled. The full model follows the relatively flat distribution 
at central rapidities shown by the data. 

In Fig. |] we compare the DPMJET-III model with and without the diquark 
breaking diagrams with data on net-proton production in central S-S collisions. 
Also here the significant dip at central rapidity in the model without the new 
diagrams is much less pronounced in the full model, however, the agreement to 
the data ||] is not perfect. 

The presence of the new baryon stopping diagrams modifies also the ex- 
trapolation of multistring models to cosmic ray energies. The energy fractions 
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Fig. 3. Net proton (p — p) rapidity distribution in p-Au collisions. The DPMJET-III 
results with and without (No Sto) the diquark breaking diagrams are compared with 
data 0. 
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Fig. 4. Net proton (p—p) rapidity distribution in central S-S collisions. The DPMJET- 
III results with and without (No Sto) the diquark breaking diagrams are compared with 
data j|. 



carried by baryons decrease against those predicted by models without the new 
diagrams. The energy fractions carried by mesons and the spectrum weighted 
moments of mesons increase as compared to models without the new diagrams. 
We present as function of the energy for two important variables the predictions 
of DPMJET-III with and without the new baryon stopping diagrams. All our 
plots are for p-p collisions, the model behaves in a rather similar way also for 
p-Air collisions. We first discuss plots, where the baryon stopping mechanism 
causes significant differences. 
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The cosmic ray spectrum-weighted moments arc defined as moments of xi ab 
distribution for the production of secondary particles i in hadron-hadron and 
hadron-nucleus collisions 

Fi(xi ab ) = xiab-j— 1 - (!) 



as follows 



fi= {xiabV 1 F i (xi ab )dxi ab . (2) 
Jo 

Here —7 ~ —1.7 is the power of the integral cosmic ray energy spectrum. The 
spectrum-weighted moments for nucleon-air collisions, as discussed in 
determine the uncorrelated fluxes of energetic particles in the atmosphere. 
We also introduce the energy fraction K t : 

Ki = / Fi(xi a b)dxi ab (3) 







As for Xi ab , the upper limit for K is 1 in h-nucleus collisions. 

In Fig.^J we present the spectrum weighted moments for charged pions in p-p 
collisions as function of the cms energy y/s. 
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Fig. 5. Spectrum weighted moments for charged pion production in p-p collisions as 
function of the cms energy ^/s. 



In Fig.^j we present again for pp collisions the energy fraction K for net 
baryons B — B (baryon minus antibaryon). The difference between K b _q and 
Kb is the energy fraction going into antibaryons Kg which is equal to the energy 
fraction carried by the baryons which are newly produced in baryon-antibaryon 
pairs. 
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Fig. 6. Laboratory energy fractions for net baryons (baryon minus antibaryon) B — B 
production in p-p collisions as function of the cms energy y/s. 



There are also observables where the difference between the two versions 
of the model are rather insignificant. Examples are the average transverse mo- 
mentum of charged hadrons as function of the energy and the average charged 
multiplicity < n c h > as function of the collision energy. 
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